Many-body calculation of the spatial extent of the wave-function of a non-magnetic 
impurity in a d-wave high-temperature superconductor using the t— J model 
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Scanning tunneling microscopy (STM) by providing images of the effects of individual zinc- 
impurities in cuprate superconductors with unprecedented atomic-resolution offers a stringent test 
to models of correlated fermions for high-temperature superconductors. Using a t-J model sup- 
plemented by Variational Monte Carlo many-body techniques, the spatial dependence of the hole 
density and of the valence bond and superconducting pairing amplitudes around the impurity are 
computed. A cross-shaped four-fold symmetric structure very similar to the observed STM obser- 
vation is found, giving strong credit to the model. 



Introduction - Cuprates superconductors can be con- 
sidered as doped two-dimensional (2D) Mott insulators 
where electronic correlations play a dominant role^'^. A 
number of exotic properties such as the pseudo-gap be- 
havior reflect the complexity of the system. In a pio- 
neering work, Anderson proposed the Resonating Valence 
Bond (RVB) Mott insulator as the relevant underlying 
parent stated from which gapless d-wave superconductiv- 
ity naturally emerges under doping. Within this scenario, 
the pseudo-gap naturally emerges as the energy scale as- 
sociated to the formation of singlet electron pairs via the 
nearest-neighbor antiferromagnetic (AF) exchange. A 
mean-field version of the RVB theory using a t-J model^ 
could also explain a number of bulk experimental obser- 
vations^. 



t, J 



t=j=o X t^J^o 



Local probes of correlated materials with atomic reso- 
lution have recently become possible thanks to Scanning 
Tunneling Microscopy (STM)^ which provided unprece- 
dented high-resolution maps of the surface of some under- 
doped cuprate superconductors?. The measured space- 
resolved doped-hole charge density in the SC regime 
of Na-CCOC and Dy-Bi2212 cuprates revealed stripe 
patterns'". This discovery naturally raises the question 
whether such inhomogeneities are induced by impurities 
or whether they are intrinsic as the bulk static charge 
and spin stripe orders detected in Nd-LSCO^ and LBCO^ 
cuprates at doping S ^ 1/8. 

Substituting a single impurity atom for a copper atom 
indeed strongly affects its surrounding region. Therefore, 
it can serve as a local fine probe, providing important 
insights about the properties of the correlated medium 
itself Imaging the effects of individual zinc impurity 
atoms on superconducting Bi2212 performed by STM^^ 
showed clear real-space modulations which can be con- 
fronted to theoretical modeling. In other words, such 
observations offer a new stringent test to models of cor- 
related fermions for high-temperature superconductors. 
It has been argued that a number of bulk properties of 
these materials can be explained within the correlated 
t-J modeli^. However, local real-space responses have 
not yet been calculated reliably since, due to the short 
superconducting coherence length, a fully many-body ap- 



FIG. 1: Schematic representation of the impurity model. The 
cross corresponds to the impurity site iq. The parameters t 
and J are set to zero on the four bonds connected to the 
impurity site. 



proach is needed. In this Letter, we have carried out such 
a program of (i) computing within a many-body numer- 
ical technique the response induced by the introduction 
of a spinless impurity site within a bulk two-dimensional 
t-J model and (ii) confronting the theoretical results to 
the experimental observations to validate or invalidate 
the model. 

The t — J Hamiltonian on a square lattice reads, 



h.c. 



Si • Sj 



(1) 



A zinc impurity is in the same 2-1- oxidation state as the 
copper ion it is substituted for, so that it does not intro- 
duce extra charge. However, in contract to the copper 
S=l/2 ion, the Zn^+ ion is in a spin-singlet state, inert 
magnetically. Hence, one can use a simple description : 
on the four bonds connected to the impurity site we set 
tij = and Jij = as shown on Fig. [TJ On all the other 
bonds, we set tij and Jij to the same values t and J, re- 
spectively. Although completely local such a "boundary" 
is expected to lead to a spatially- extended perturbation 
strongly affecting GS properties. 

Description of method - Variational fully-projected 
fermionic wave- funct ions ^'^ (WF) are known to incor- 
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porate very satisfactorily correlation effects of the t- 
J model. We have extended them to finite periodic 
L X L clusters containing a single impurity (in practice, 
L = 8 and L = 16). A Variational Monte Carlo (VMC) 
schemei^ is used to (i) realize a thorough optimization 
over the variational parameterai^ (see below) and (ii) 
to calculate the ground-state (GS) physical observables. 
The presence of the impurity on site iq modifies the 
Hilbert space, i.e. cj^crl^) = ^ and Ciga\'^) — 0, where l^*) 
is the GS of the system. In our Monte Carlo variational 
scheme an "impurity projector" Vig = {l — nig^){l — nigi) 
is inserted, and the impurity variational wavefunction 
is defined as |^vmc) — 'Pg'^io\D), where Vg is the 
usual Gutzwiller projector enforcing the constraint of no- 
double occupancy on the remaining — \ sites. 

Motivated by the success of the RVB theory to explain 
bulk properties^ the mean-field determinant \D) is chosen 
to be the ground-state of a mean-field Hamiltonian of 
standard BCS-type, 

+ I](A^,4Tc]i+^•C■)+MX"^a, (2) 

defined on all of the L x L sites, including the io site 
(always occupied by a hole). We optimize all differ- 
ent non-equivalent bonds around the impurity, starting 
from an initial guess respecting or not the square lattice 
symmetry around iq. Since in principle the Hamiltonian 
symmetry around io (see Figure [1]) could be spon- 
taneously broken, we have performed a number of pre- 
liminary tests on small 8x8 lattices, choosing the initial 
RVB bonds pattern with lower symmetries like e.g. C4 
or C21, symmetries (the later allowing the formation of a 
domain wall). We have found that the full C^^ symme- 
try is systematically restored at the variational minimum. 
Therefore, to reduce the number of variational parame- 
ters and gain accuracy, the C^y symmetry has been en- 
forced on our largest 16 x 16 cluster. As expected, all 
optimized WFs are found to show opposite signs of Aij 
on any site-sharing vertical and horizontal bonds, hence 
reflecting the expected orbital d-wave character of the su- 
perconducting order. Lastly, we note that allowing finite 
values of the parameters Aij and Xij on the four bonds 
connected to the impurity is also important to gain en- 
ergy as shown on Table I. The lowest-energy state (II) is 
obtained for a full optimization of the Ay and Xij 
rameters over all bonds. Typically, Xij has a significant 
magnitude on the four bonds connected to the impurity 
site. Moreover, for decreasing doping, a sizable ampli- 
tude of Aij also appears on the later bonds. 

Results on 16 x 16 clusters - We now turn to the VMC 
calculations on the 16x 16 cluster with periodic-boundary 
conditions, assuming a physical value of t/ J = 3. Here, 
we consider a physical " core" 8x8 region centered around 
the impurity (i.e. of the same size as our previous 
small cluster), where we impose a C^y symmetry around 



TABLE L Variational energy per site (in units of t) for differ- 
ent projected WFs for the t — J model at doping 1/8 {Nh ~ 8), 
for t/J = 3 and a 8 x 8 cluster with an impurity. (I) 1D(8 x 8)) 
is optimized fixing Aij — 0, Xij ~ around the impurity and 
fiig — 0. Everywhere else Aij optimized and Xij = Ij (H) 
\D{8 X 8)) is fully optimized with all possible Aij and Xij 
(including the impurity bonds) . The (total) energy difference 
between (I) and (II) is ~ 0.3t per impurity. 

T?>: (I) (11) 

EvMc[t]: -0.41493(5) -0.41968(5) 

the impurity site. Outside, we assume a uniform d- 
wave superconducting background (bg) whose parame- 
ters Xi,i+x = Xi,i+y = Xhg and Ai^i+S: = —Ai^i+y = Abg 

are optimized simultaneously. This enables to reduce sig- 
nificantly the boundary effects and is justified since the 
spatial extension of the effect of the impurity rarely ex- 
ceeds the assumed size of the core. 

The spacial distribution of the local hole density 
(cio-c|g.), the bond hole kinetic amplitudes Kij ~ 
{{cja-cl^ + h.c.)) and the magnetic VB amplitudes Sij — 
(Si • Sj) are shown in Fig. ^a.) and in Fig. El^a), re- 
spectively, for doping 1/8. Here, and throughout the pa- 
per, we only show the 6x6 central region exhibiting the 
largest modulations. It turns out that the variational 
parameters Aij are suppressed on the four bonds con- 
nected to the impurity. This is compensated by an in- 
crease in Aij and hence of Sij on the neighboring bonds, 
forming a cross-like structure (see thick blue bonds of 
Fig-EIa)). Due to similarities with work done in a some- 
what different context^i, we shall refer to these bonds 
as the four "dimer bonds". These bonds are character- 
ized by a simultaneous hole deficiency and a large gain in 
the magnetic energy (which can reach more than 40%), 
hence signaling a tendency towards singlet crystallization 
around the impurity. The distribution of Kij in Fig. [2] 
shows also a remarkably strong modulation around the 
impurity. 

Superconducting properties of RVB states are char- 
acterized by the singlet-pair correlations at distance r, 
(^'vMc|A^+rAs|1'vMc)/(^'vMc|*VMc), where the oper- 
ator Aj, — c\ . I ~ c V \ , -\ t creates a sin- 
glet pair of electrons on the bond between locations s 
and s -I- a on the lattice, a being the unit vector that 
specifies the bond direction (along x or y). On the 8x8 
cluster, we have computed pairing correlations between 
separate bonds for increasing bond separation. However, 
at the largest distance available on this cluster, the cor- 
relations have not completely reached saturation. To get 
a better estimation of the superconducting order param- 
eter we have considered the 16 x 16 cluster and computed 
the pairing amplitudes Ay for all bonds (?, j) within the 
"core" region, 

(AsAb.) 

Ai(s)j(s+a) = , ° = , (3) 
V(^bgAbg) 
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FIG. 2: (Color online.) VMC results for the GS on-site hole 
densities (circles) and kinetic bond amplitudes (colored seg- 
ments) obtained on a 16 x 16 cluster. Only the central re- 
gion around the impurity is shown. Diameters of circles and 
widths of segments scale with the absolute value of the rela- 
tive differences w.r.t. the impurity-free homogenous state at 
the same 5avc hole density (whose reference values are esti- 
mated by interpolating pure clusters with available flanking 
hole densities) . Higher (lower) hole densities and bond magni- 
tudes w.r.t. the homogeneous case are shown by open (filled) 
circles and blue (green) bonds respectively. For complete- 
ness, we also show on the plot the (bare) numerical values of 
the non-equivalent sites/bonds, (a) and (b) corresponds to 
32 ((5avo ^ 0.1255) and 20 ((5ave ^ 0.0784) doped holes giving 
rise, for an homogeneous background, to E^°'^°^ /t = —0.1487 
and E^°^°^ /t = —0.0941 per bond, respectively. 



where Abg is a pair operator on the most remote bond 
in the homogeneous background. As shown in Fig. ID^a), 
pairing is enhanced on the dimer bonds, and is depleted 
around the impurity, where holes are less present. 

The 16x16 cluster also allows to reduce the doping 
content, e.g. to 20 holes, going further into the under- 
doped region. Interestingly, the hole distribution around 
the defect is very sensitive to the doping ratio. Indeed, 
for doping around 12.5% (32 holes) we found that holes 
are slightly repelled from the bonds around the impurity. 
In contrast, for 7.8% doping, holes tend to concentrate 
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FIG. 3: (Color online.) VMC results for the GS magnetic 
bond amplitudes obtained on a 16 x 16 cluster. Same conven- 
tions and parameters as in Fig[2] (a) and (b) corresponds to 
32 (5avo =i 0.1255) and 20 (5avc ==; 0.0784) doped holes giving 
rise for an homogeneous background to E^^°^/J = —0.074 
and E^^°<^/J = -0.084 per bond, respectively. 



more around the impurity site as shown in Fig.[2jb). The 
variational pairing becomes stronger on the impurity 
bonds suggesting the formation of a "hole pair" with the 
impurity empty site. The corresponding real space mod- 
ulations of Sij and Ay are shown in Figs.[3]^b) and[3I^b). 

Discussions - Let us now compare our findings to prior 
theoretical approaches. The first investigation of a sin- 
gle impurity immersed in a correlated host has been per- 
formed using Lanczos exact diagonalization of small clus- 
ters. A calculation of the local density of states^^ re- 
vealed bound-states (of different orbital symmetries) in 
which a mobile hole is trapped by the induced impurity 
potential. Here, a unique mobile hole was assumed in the 
cluster, hence preventing real bulk pairing and giving rise 
to a very small doping ~ 5% in the surrounding region. 
Although our VMC calculations are done in a different 
physical range (and on much larger clusters), we find, for 
decreasing doping, the emergence of excess hole density 
around the impurity, which possibly could be consistent 
with a bound-state formation when 5 — > 0. 

Metlitski and Sachdev^^ have introduced a theory of 
valence bond solid (VBS) correlations near a single impu- 
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FIG. 4: (Color online.) VMC results for the pairing bond 
amplitudes Kij obtained on a 16 x 16 cluster. Same conven- 
tions and parameters as in Fig[2] (a) and (b) corresponds to 
32 (5avc ^ 0.1255) and 20 (5ave ^ 0.0784) doped holes giving 
rise, for an homogeneous background, to Abg — 0.0787 and 
Abg ^ 0.0626, respectively. 



rity in a square lattice antiferromagnet. When the system 
is close to a quantum transition from a magnetically or- 
dered Neel state to a spin-gap state with long-range VBS 
order, a missing spin gives rise to a VBS pinwheel (or 
"vortex") around the impurity^^. To compare w\t\v these 
predictions we have computed the VBS order parameter 
of Eq.(5) in Ref. [TtI . However, despite many similari- 
ties (e.g. crystallization of dimer bonds in the vicinity 
of the impurity) , the vortex structure is not recognizable 
in our simulation. We hypothesize that (i) our system is 
probably not close enough to the critical point assumed 
in Ref. ll6lfT7l and/or (ii) the VBS region develops differ- 
ently in a d-wave RVB than in an AF background. 

Lastly, our results are compared to the experimental 
STM observations around a Zn impurity in a Bi2212 
cuprate superconductor shown in Ref. illi. First, we point 
out that GS properties have been calculated here while 
Ref. [Til reports spectral properties. However, since equal- 
time and frequency-dependent quantities are related via 
a simple frequency integration up to a physical cutoff 
(see e.g. Ref. for derivation of the local hole-charge 
distribution from space-resolved tunneling spectra) both 
sets of data should reveal similarities. Indeed, like in the 
experiments, the patterns we found show clearly a cross- 
shaped symmetric structure providing evidence that (i) 
the theoretical modeling of the impurity as an empty site 
and (ii) the use of the strongly correlated t-J model to 
describe the bulk high-Tc superconductor are realistic. 
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